Probing dark energy with gamma-ray bursts 
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We propose a new method to use gamma-ray bursts (GRBs) as an alternative probe of the dark 
energy. By calibrating luminosity- variability and luminosity-lag time relations at low redshift where 
distance-redshift relations have been already determined from type la supernova (SNIa), GRBs at 
high redshift can be used as a distance indicator. We investigate the potential impact of future GRB 
data on determining the current matter density Q m and the dark energy equation of state w which 
is assumed to be constant. We show that a combined analysis of a data set expected from the Swift 
and the current SNIa data results in excellent determination of both Q m and w. 



Accelerating universe, which is strongly suggested by 
the distance measurements of type la supernova (SNIa) 
[1-3] , is one of the most interesting and challenging prob- 
lems in modern cosmology. Recent observation of the 
cosmic microwave background (CMB) anisotropy [4] re- 
vealed that the universe is almost flat and confirmed the 
need for the energy component which accelerates our uni- 
verse. 

The cosmological constant, or "dark energy" in gen- 
eral, is often introduced to account for this cosmic accel- 
eration. Dark energy is the unknown energy component 
which has effective negative pressure p = wp (w < 0). 
The case with w — — 1 corresponds to the cosmological 
constant. Since the value of w, which is a function of 
redshift z in general, is closely related to the nature of 
dark energy, it is of great importance to determine the 
value of w observationally. Until now, a number of meth- 
ods are proposed to probe the dark energy equation of 
state [5-13]. However, currently w has not been deter- 
mined quite well even if one assumes that w is constant 
[14,15]. In particular, the current distance measurements 
of SNIa poorly determines the value of us because of the 
degeneracy between w and the dark-energy density [1-3]. 
Moreover, the use of CMB anisotropy, which is known to 
be one of the most powerful and clean probe of cosmolog- 
ical parameters, introduces an additional degeneracy be- 
tween the Hubble constant and dark-energy parameters 
[15]. That is one of the reason why future SNIa search, 
such as the proposed SNAP satellite [16], has been vig- 
orously studied as a probe of the dark energy [17-19]. In 
this Letter, we propose another method using gamma-ray 
bursts (GRBs) as a distance indicator. 

Recent studies have suggested the possibilities that 
GRBs can be standard candles through the relations be- 
tween absolute "isotropic-equivalent" luminosity (L) and 
other observed quantities. Two of them, which we use in 
this Letter, are variability (V) [20,21] and spectral lag 
(Ti a g) [22,23]. The variability is a specific measure of the 
"spikiness" of the light curve. The spectral lag is the time 
between peaks as recorded at high and low photon ener- 



gies. These two quantities are related to the luminosity 
as L cx V av and L oc , respectively. Furthermore, a 
l//riag relation, which must exist if the L/V and L/r\ ag 
relations are true, has been confirmed with an indepen- 
dent sample of 112 BATSE bursts [24]. There are also 
some theoretical explanations for these relations [25,26]. 
Once these relations are established by observations, then 
one can use these relations to derive the absolute lumi- 
nosity of GRBs and to estimate the luminosity distance 
to the GRBs [27]. 

Up to now, only ~ 10 GRBs have the enough informa- 
tion (z, V, Tiag) to establish the above two relations. How- 
ever the Swift satellite [28], planned for launch in 2003, 
is expected to detect ~ 1000 GRBs, about half of which 
will be observed with known redshifts, during 3 years of 
observations. This huge number of GRBs allows us to 
probe the expansion history of the universe more deeply 
than SNIa [27] . In this Letter, we simulate the data which 
is expected to be obtained by the Swift and investigate 
the ability to probe the dark energy independently with 
other methods such as CMB anisotropy. Throughout the 
Letter, we assume a fiat universe fi m + fix = 1, where 
f2 m is the present matter density and fix is the present 
dark-energy density. 

Here we must be careful to avoid circular logic. The 
original studies of the L/V and £/ T iag relations assumed 
a particular set of cosmological parameters (the Hubble 
constant Hq, w, fl m and the flatness of the universe) to 
derive the relations. Without any assumptions about the 
cosmological parameters, we cannot determine L/V and 
L/n ag relations which are used to obtain the absolute 
luminosity of GRBs. Thus we assume that a form of 
the luminosity distance as a function of the redshift at 
2 < 1 has already been well determined by the current 
SNIa observations. This is approximately that of the A- 
dominated universe (w — — l,flx — 0.7 and fl m = 0.3). 
The errors associated with the form of the luminosity dis- 
tance at z < 1 is reflected in the systematic errors in the 
derived absolute luminosities of the low-z [z < 1) GRBs. 
We calibrate these two relations by using only low-z 
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GRBs, and obtain the absolute luminosities of high-z 
(z > 1) GRBs using these relations. After that, the pro- 
cedure is the same as that using SNIa: fitting parameters 
such as w and fi m so as to reproduce the magnitude as a 
function of the redshift. The apparent magnitude m and 
redshift z are related by the luminosity distance as: 



(1) 
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x [O m (l + z'f + (1 - fl m )(l + z'f +3w } ~ h (2) 
is the luminosity distance in units of Mpc. 



m = M + 25 + 51ogd L (z), 
where M is an absolute magnitude and 

d L { z ) = (1 t z)c / ,/-.' 




FIG. 1. Redshift distribution of the simulated GRBs. The 
flux limit of the Swift is applied. 

For concreteness, we simulate a future data set similar 
to one expected from the Swift. The sample comprises 
500 GRBs, which is consistent with the expected num- 
ber of GRBs with measured redshifts. Their redshifts 
and luminosities are distributed according to the GRB 
rate history and the luminosity function based on the 
model of the star formation rate 2 (SF2) in [29], which is 
roughly constant at z > 2. This model is consistent with 
the current observations of UV luminosity density of the 
galaxy population. The energy spectrum of the GRB is 
assumed to be a single power-law with index —2.25. The 
observed luminosities are calculated assuming a back- 
ground cosmology with Ho — 65kms _1 Mpc~ 1 , Sl m = 
0.3, fix = 0.7 and w = -1. The flux limit of the Swift, 
f > 0.04 photons cm _2 s _1 , is applied to check whether 
a produced GRB will be observed or not. The errors in 
redshifts are tentatively neglected because of the lack of 
information on the redshift measurement by Swift, while 
we will show later that inclusion of the errors in redshifts 
does not change our result so much. We assume the log- 
normal distribution of the luminosity errors of 0.02 in 
logarithmic units of base 10. Other errors below are also 
assumed to be the log-normal distribution. Fig. 1 shows 



the redshift distribution of the simulated GRBs. The av- 
erage redshift of these GRBs is about 3. Although very 
high-z GRBs might be subject to some possible biases, 
they are practically not influential on the analyses below, 
because the number of such high-z GRBs is not so large. 
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FIG. 2. Absolute luminosity-variability relation for the 
simulated data. Points with error bars show low-z [z < 1) 
GRBs, which are used to derive L/V relation. Points with- 
out errorbars show high-z (z > 1) GRBs. Note that the high-z 
GRBs also have the same errors as low-z GRBs. Solid line 
shows the assumed L/V relation (3). 
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FIG. 3. Same as Fig. 2, but the absolute luminos- 
ity-spectral lag relation is plotted. Solid line shows the as- 
sumed L/r\ ag relation (4). 

Next the observed variability and spectral lag are given 
to each GRB according to the following L/V and L/r\ ag 
relations, 



L = 3 x 10 56 (V(l + z)) erg/s, 

.so (T\ ag {i + zy _ ! 
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with the intrinsic scatter as well as the measurement un- 
certainties. For both the L/V and L/r\ ag relations, we 
assume both the intrinsic scatter and measurement error 
to be 0.1 in logarithmic units. The indices and propor- 
tional constants of these relations are different among 
various studies. The values we use are typical of them 
[20-23]. The measurement errors are almost the same 
as those of the current observations. The simulated data 
are plotted in Fig. 2 and 3 with the assumed relations (3) 
and (4), respectively. Although the L/r\ ag relation might 
steepen to L oc at T\ ag > 0.1s [23], our assumption 
of single power-law for the L/r\ ag relation is enough be- 
cause GRB events with such large ri ag will be very rare 
for the luminosity function adopted here. In fact the in- 
dices that we adopted need not be the true values since 
they are to be determined from the Swift observation. 
Thus what we require are that there are some one-to-one 
relations between the observable quantities (V and Ti ag ) 
and the absolute luminosity L and that the relations are 
independent on z for < z < 3. Although evolution of 
the relations is not likely as is stated in [27], it will be 
tested if the luminosity distance up to z < 1.7 is obtained 
from the SNAP observation. 

Then the absolute luminosities of the low-z GRBs are 
calculated by the luminosity distance-redshift relation 
which has been determined by the current SNIa obser- 
vations. Its uncertainties lead to the uncertainties of the 
absolute luminosities of ~ 0.2 in logarithmic units. Then 
the L/V and L/r\ ag relations are derived by the low-z 
GRBs. As a result, we recover the the L/V and L/r\ ag 
relations of the form L oc V 213 and L oc from 
the simulated low-z GRBs. Two independent absolute 
luminosities of high-z GRBs are obtained by these rela- 
tions, and are combined as weighted averages to produce 
combined absolute luminosities. Fig. 4 shows the binned 
magnitude difference Am between several representative 
cosmological models and the fiducial model (£! m = 0.3, 
w = — 1), along with the simulated data. 
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FIG. 4. Magnitude difference Am between several cosmo- 
logical models and the fiducial model (fl m = 0.3, w = —1), 
along with the simulated binned data. Three solid curves 
are for (Cl m , w) = (0.2, -1), (0.3, -1) and (0.4, -1), and two 
dashed curves are for (Q m ,w) = (0.3, —1.5) and (0.3, —0.5). 



Finally we perform \ 2 constraints to see to what extent 
we can constrain dark energy properties from the obser- 
vation of GRBs. Here we assume the dark energy to have 
a time- independent equation of state w = p/ p = const. 
We do not put a prior w > — 1 both because it may bias 
our final result [14] and because several theoretical mod- 
els which predict w < — 1 do exist [30,31]. Further a 
flat universe with Hq = 65kms^ 1 Mpc~ 1 is assumed. We 
calculate x 2 as 
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where Agrb is the number of GRBs, to; is the observed 
apparent magnitude, M; is the absolute magnitude es- 
timated from L/V and L/r\ ag relations, and <5to; is the 
error of the magnitude. SNIa data fit and the combined 
fit of GRB and SNIa data are also performed. For the 
SNIa data, we use the data of the SCP (primary data 
set C) [2], the High-Z Search Team [3] and Calan-Tololo 
survey [32] for the SNIa analysis. Fig. 5 shows the best- 
fit ler and 3cj confidence regions in the Q m — w plane 
by the simulated GRB data, the current SNIa data and 
the combined analysis. Since average redshifts are signif- 
icantly different between SNIa (~ 0.5) and GRB (~ 3), 
the combined fit results in the excellent determination of 
both fl m and w. Our simulation demonstrates that Q m 
and w can be determined within ~ 20% accuracy from 
the combined analysis of only two types of observation;, 
SNIa and GRBs. Note that constraints from both SNIa 
and GRBs are quite insensitive to other (cosmological) 
parameters, such as the Hubble constant, the spectral in- 
dex of primordial power spectrum, and the normalization 
of rms mass fluctuations, which are also very important 
in using other methods, e.g., CMB anisotropy or cluster 
abundance. 
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FIG. 5. Best fit la and 3<r confidence regions in the Q m — w 
plane by the simulated GRB data (dotted), the current SNIa 
data (dashed) and the combined analysis (solid). The cross 
shows the assumed model (fl m ,w) = (0.3, —1). A flat universe 
is assumed. 
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Besides the fiducial simulation, we perform two addi- 
tional pessimistic simulations: one with 300 GRBs in- 
stead of fiducial 500 GRBs, and another with redshift 
error of 0.2 in logarithmic units which is neglected in 
Fig. 5. The best fit la confidence regions for these two 
pessimistic simulations as well as the fiducial simulation 
are shown in Fig. 6. Although the two pessimistic sim- 
ulations result in larger la regions, tt m and w are still 
determined well. Therefore, we conclude that search for 
distant GRBs can be a promising way to probe the dark 
energy density and dark energy equation of state sepa- 
rately. 
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FIG. 6. Best fit la confidence regions of the combined anal- 
ysis for two additional pessimistic simulations as well as the 
fiducial simulation. Solid, dotted and dashed curves cor- 
respond to the result of the fiducial simulation, one with 
300 GRBs, and one with redshift error of 0.2 in logarith- 
mic units, respectively. The cross shows the assumed model 

(n m , w ) = (o.3,-i). 

In summary, we have investigated a potential impact of 
future GRB searches on probing the cosmological param- 
eters such as the matter density fraction of the universe 
fl m and the equation of state of dark energy w. Our re- 
sult is that a combined analysis of a data set expected 
from the Swift and the current SNIa data results in the 
excellent determination of these parameters. 
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